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A simple technique is described that provides an improved nearside-farside (NF) method for elastic 
scattering amplitudes. The technique, involving the novel resummation of a Legendre partial wave 
series, reduces the importance of unphysical contributions to NF subamplitudes, which can arise 
in more conventional NF decompositions. Detailed applications are made to a strong absorption 
model and to a ^^O + ^^O optical potential at Siab = 145 MeV. We also discuss ^"^O + ^^O at £iab 
= 480, 704, 1120 MeV, and a + ^^C, a + "'"Ca, both at Biab = 1370 MeV. 

PACS numbers: 24.10.Ht, 25.70.Bc, 34.50.-s, 03.65. Sq 



I. INTRODUCTION 

In heavy-ion, atomic and molecular collisions, an elas- 
tic differential cross section, a{6), where 9 is the scat- 
tering angle, is often characterized by a complicated in- 
terference pattern. This complicated structure makes it 
difficult to understand the physical phenomena involved 
in the scattering process, as well as the links between 
(j{9) and the properties of the model that describes the 
phenomenon. 

In some cases, semiclassical methods Q explain the 
scattering pattern as the interference between simpler, 
and slowly varying, subamplitudes. If we ignore the com- 
plication that, in some angular regions, uniform asymp- 
totic techniques are necessary, then the semiclassical sub- 
amplitudes arise mathematically from saddle points or 
poles, which account physically for contributions from 
reflected, refracted or generalized diffracted semiclassical 
trajectories The subamplitudes can be conveniently 
grouped into two types: those arising from semiclassical 
trajectories which initially move in the same half plane as 
the detector (N or nearside trajectories) and those from 
the opposite half plane (F or farside trajectories). 

Semiclassical methods are not always simple to apply 
and sometimes they have a limited range of applicability. 
Their limitations are determined by the range of validity 
of the (presently known) asymptotic techniques that are 
used to approximate the original quantum mechanical 
problem. 

In order to overcome these difficulties, it is common 
practice to apply to the elastic scattering amplitude, 
f{9), a NF method that was proposed by Fuller ||] more 
than 25 years ago. The Fuller NF method has the merit 
of being simple and, although inspired by the semiclas- 
sical theories, it uses only scattering matrix elements Si 
calculated (or directly parametrized) by exact quantum 



'Electronic address: 
t Electronic address: 



R.aimondo.Anni@le.infn.it 



J.N.L.Connor@Manchester.ac.uk 



mechanics. The NF subamplitudes are obtained by ex- 
act summation of NF partial wave series (PWS), thereby 
bypassing problems associated with the applicability and 
validity of semiclassical techniques, e.g. using approxi- 
mate Si, replacing the PWS by integrals, using stationary 
phase integration, etc. 

The Fuller method is based on a splitting of the Legen- 
dre polynomials, P/(cos0), in the PWS of f{9), into trav- 
eling angular wave components, with the traveling angu- 
lar waves consistent with detailed semiclassical analyses 
of scattering from impenetrable and transparent spheres 

0- 

The Fuller method was || and continues to be 
widely used; indeed the ISI Web of Science reports about 
140 citations since 1981 to Fuller's seminal work (for more 
recent examples see [|[ 0, ^, ^ 11 and references 
therein). The success of the method depends, apart from 
its simplicity, on its remarkable capability of physically 
explaining complicated interference patterns in cross sec- 
tions as arising from the interference between NF subam- 
plitudes having simple properties. In particular, the NF 
cross sections are often less structured and more slowly 
varying with 9 than is the full cross section. Even though 
no semiclassical technique has been used, these NF sub- 
amplitudes can often be given a physical interpretation 
(in analogy with results from semiclassical methods) as 
contributions from simple scattering mechanisms, which 
then allows a good understanding of the angular scatter- 
ing. 

In the light of these unquestionable successes, it is de- 
sirable to extend the validity of the Fuller approach to 
cases where the original Fuller method is not (physically) 
satisfactory, for example, it may produce oscillatory and 
rapidly varying NF cross sections, when the full cross 
section is monotonic and slowly varying with 9. Exam- 
ples of these shortcomings have been known for a long 
time. One classic example is pure Coulomb scattering. 
For repulsive Coulomb potentials only a N contribution is 
expected semiclassically ([|l|, p. 56), whereas the Fuller 
NF method yields also a F contribution As a re- 
sult, the NF cross sections are less simple than the full 
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one. In this case, the unsatisfactorily effects are, however, 
confined to a restricted backward angular region (@, p. 
1564). Another more striking example is observed in the 
scattering by a uniformly charged sphere (Q, p. 154, 
Fig. 26). In this case, the ratio of the full cross section 
to the Rutherford one decreases monotonically into the 
shadow of the Coulomb rainbow. In contrast, the N ra- 
tio closely follows the full ratio up to 6* ~ 40° when it 
becomes approximately constant (i.e. independent of 9), 
being approximately equal to the F ratio. 

A similar effect is also observed in the angular dis- 
tributions for a strong absorption model (SAM) with a 
two parameter (A and A) symmetric S matrix element 
and Fermi- like form factors For a fixed value of 

the cut-off parameter A and for a sufRciently large value 
of the diffuseness parameter A, the Fuller NF cross sec- 
tions show an almost exponential decline up to a certain 
6 (which decreases with increasing A). At larger angles, 
the NF cross sections are greater then the full cross sec- 
tion, which continues its oscillatory exponential decline. 

Similar striking effects appear at high energies and for 
large scattering angles in the NF cross sections of a parti- 
cles and light heavy-ions scattered by nuclei using optical 
potentials. Less striking, but still disturbing, effects are 
also observed, at lower energies, in some N cross sec- 
tions for the optical potentials used to fit recent data of 
light heavy- ion scattering [|[ |[ H, |l^. Typical N cross 
sections rapidly decrease from 0° and from 180°. The 
two branches meet in the crossing region where an inter- 
ference pattern appears, with strong oscillations over an 
extended range of angles. 

In this paper, we show how some of these shortcom- 
ings can be removed using a new NF method , based 
on an improved modified resummation of the PWS. The 
new method is a development of the Hatchell idea of 
incorporating the Yennie, Ravenhall and Wilson (YRW) 
[Q resummation technique into the NF formalism. The 
limitations of the NF Hatchell resummation technique 
have been discussed |15, 1^, |l^ and a modified NF YRW 
resummation procedure, depending on two parameters a 
and /3, was proposed ||l^, |l9| to bypass the difficulties 
with the original NF YRW approach. The possibility of 
further improving the modified YRW resummation proce- 
dure, using different resummation parameters ai,a2, ■ ■ ■ 
and /3i,/?2, . . . , together with a rule to fix the value of 
these parameters [|l3|, is discussed in the present work. 

For all three NF methods. Fuller, Hatchell and ours, 
the starting point is the quantum mechanical PWS for 
the full scattering amphtude f{9), 



(=0 



(1) 



In the following, we will write x = cos 9. For future 
convenience, we recall that the PWS (|l|), considered as 
a distribution, converges to its exact value for 9^0, 
upon dropping the 1 in the term S*; — 1 on the r.h.s of 
(^). The omitted amplitude is proportional to the Dirac 
delta function S{1 — x) (e.g. see p. 52). We also 
recall that the PWS (|l|), considered as a distribution, is 
convergent if Si is asymptotically Coulombic |Q . 

In Sec. ||, we briefly outline the original Fuller NF 
method, and show that unsatisfactorily results are ob- 
tained for seven collision systems, which are different 
from those considered in [flSl. In Sec. IH, we discuss 



a modification of the Fuller NF technique proposed in 
©0 



Refs. 
method 
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|19[ and present an improved modified 
Our new method is very effective in cleaning 
unphysical contributions from NF cross sections for the 
seven examples where the usual Fuller technique gives 



usatisfactorily results. Our conclusions are in Sec. IV 



II. LIMITATIONS OF THE FULLER NF 
METHOD 

A. Introduction 

The Fuller NF decomposition is realized by splitting 
in (|^) Pi{x), considered as a standing angular wave, into 
traveling angular wave components 



Pi{x)^q[-\x) + q\+\x), 
where (for a; ^ ±1) 

Q\^\x)^hPiix)±-Qiix)], 

Z TT 



(3) 



(4) 



with Qi{x) the Legendre function of the second kind of 
degree I. 

Inserting (j|) into (j^) splits f{9) into the sum of two 
subamplitudes 



with 



f{9) = f^-\9)+f^+\9), 



(5) 



(6) 



(=0 



Note that, by construction, the decomposition (||) is ex- 
act. Also we will obtain an exact decomposition, by using 
in place of the Q^i^\x) in (^), any pair of functions whose 

sum is Pi{x). The property of the q')^\x) that makes 
the splitting (0) important, is the asymptotic result 



where k is the wavenumber, Pi{cos9) is the Legendre 
polynomial of degree I and a; is given in terms of the 
scattering matrix element Si by: 



ai = {2l + l){Si-l). 



(2) 



gp^(x) 



for IsmO ^ 1, where A = / -|- ^. In particular, (^ allows 
(— ) to be identified with N scattering and (+) with F 
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scattering (Q, p. 121). In the semiclassical theory, the 
spHtting of Px-^ix) into the sum of Q'^^iix), or the re- 
lated sphtting obtained from the asymptotic expansions 
of these functions plays a crucial role in deriving 
the semiclassical subamplitudes. In particular, the NF 
semiclassical subamplitudes arise from terms originally 
containing Q''^i{x), or their asymptotic expansions (|^). 
These facts raise the hope that the direct calculation of 
the /'-^-'(fi') from their PWS representation in (^) will 
separate the NF contributions to f{9), thereby avoiding 
problems connected with the applicability or validity of 
the semiclassical theory. 

In order to make this hope mathematically rigorous, 
one should prove that it is possible to perform on the 
PWS, written in terms of the Q[^\x), the same manip- 
ulations that are used in deriving the complete semiclas- 
sical decomposition of f{9). These manipulations are 
essentially path deformations in A of the integrals into 
which (|l|) can be transformed, using either the Poisson 
sum formula ([Q, p. 45) or the Watson transformation 
([Q,p. 49). The consequences of these path deformations 
depend on the properties of the terms in the PWS when 
they are continued to real or complex values of A from 
the initial physical half integer A values. The splitting 
of Piix) into the qI'^\x) modifies these properties and 
can cause the appearance of unphysical contributions in 
the f^^^{9), which cancel out in f{6) (these contribu- 
tions are not expected to be present in the semiclassical 
subamplitudes) . 

In spite of these possible limitations, extensive experi- 
ence with the Fuller NF method has demonstrated that 
the method is usually reliable, in the sense that it often 
decomposes f{9) into simpler NF subamplitudes, appar- 
ently free from unphysical contributions arising from the 
above mathematical difficulties. However for a few exam- 
ples, some of which were mentioned in Sec. I, the Fuller 
NF subamplitudes can be directly compared with the cor- 
responding exact analytical, or semiclassical, results and 
disagreements are observed. 

Fortunately, the Fuller NF subamplitudes contain in- 
formation that allows one to recognize the unphysical na- 
ture of the undesired contributions. Suppose Z'^^^), or 
/^~-'(0), contains a single semiclassical contribution from 
a stationary phase point at X{9). Then the derivative 
with respect to 9 of the phase of /^^H^), or f^~^{9), is 
equal to X{9), or —X{9), respectively ([Q, p. 57). Follow- 
ing Fuller we will call this derivative the Local Angular 
Momentum (LAM) for the F (or N) subamplitude; it de- 
pends on 9. Usually, only for a generalized diffracted tra- 
jectory, arising from a simple pole, is the LAM expected 
to be constant, being equal to the angular momentum of 
the incoming particle responsible for the diffraction. 

In the semiclassical regime, this constant value is ex- 
pected to be large. Because of this, if we observe that, 
in a certain 9 range, LAM w 0, this can be considered a 
warning for the unphysical nature of the N or F subam- 
plitudes in that range of 9. This occurs for the LAM of 



the Fuller Coulomb F subamplitude, and for the NF sub- 
amplitudes of the SAM in the angular region where the 
NF cross sections contain unphysical contributions [ p^ . 
In both cases, this decoupling of 9 from LAM, together 
with the fact that the full cross section is simpler then 
the NF ones, suggests the unphysical nature of the NF 
subamplitudes . 

We show below for seven collision systems, different 
from those considered in , how unphysical NF contri- 
butions manifest themselves. The seven collision systems 
are; (a) a simple SAM model, (b) ^^O -t- ^^O at Ei^b = 
145 MeV, using the WS2 optical potential of Ref. 0, 
and (c) more briefly, ^^O + ^^O at £'iab = 480, 704, 1120 
MeV, and a + ^^C, a + '^^Ca, both at Ei^h = 1370 MeV. 



B. Strong absorption model for elastic scattering 

The first example is a simple SAM in which the Si is 
directly parametrized by 



^'-^^^^"l + exi(A^) + l-Hexp(AJA)' 



(8) 



with A = Z i, A = 10.0 and A = 1.8, 

The SAM (with modifications for the Coulomb inter- 
action, and in a slightly different form from (j^) ) was 
widely used in early studies of heavy-ion elastic scattering 
pl| , p2| . At the present time, the SAM is not so popular, 
using either simple forms such as (H) or more sophisti- 
cated functions. It has been found that the character- 
istics of heavy-ion angular distributions, measured over 
wide angular ranges, are not accounted for by the sim- 
pler SAM models; instead the angular scattering is more 
easily described using optical potentials, rather then at- 
tempting complicated extensions of the SAM. 

In spite of this, the SAM in its simple form (^) con- 
tinues to be of interest, since it allows important tests 
of NF decompositions (l2l |6[ 0. This is because 
the PWS for the SAM ( pi) can be evaluated easily by 
saddle point techniques or more simply, using the 
Watson transformation and elementary complex integra- 
tion. Both methods allow a simple, mathematically cor- 
rect, identification of the NF subamplitudes. For A ^ 1, 
exp(— 27rA) <C 1 and Asin6' 3> 1, it is found that the 
NF subamplitudes are oc exp(--7rA0 =F i\9)/\/s\n9, to a 
good degree of approximation (Ref. [Q, Eq. (3.5) ). The 
NF cross sections, multiplied by sin^ are equal and have 
an exponential slope, whereas the phase derivatives of 
the NF subamplitudes are expected to be equal to =fA, 
respectively. 

The results obtained by applying the Fuller NF method 
to the SAM with parameters A = 10.0 and A = 1.8 
are shown in Fig. 1. In the lower panel, we show 
a log plot of the dimensionless quantity 4k^a{9) sin 9 
versus 9 since the corresponding NF quantities are ex- 
pected to have an exponential slope. This is additionally 
shown in Fig.l by the thin dot-dash line which repre- 
sents log]^Q[exp(— 27rA6')]. Furthermore, because the Si 
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are real, f{9) has a constant phase (and its phase deriva- 
tive is of no interest), while the f^^^d) have identical 
moduli but opposite phases. Thus we need only show 
the N (or F) LAM and similarly for the cross sections. 
In Fig. 1, the N and F quantities are shown by thick 
continuous and dashed curves, respectively. 

In a systematic notation explained in Sect. Ill, the 
results obtained from the Fuller NF method, which sub- 
stitutes (^ into (0), are indicated by R = 0. The thin 
curve, in the lower panel of Fig. 1, shows the full cross 
section. Figure 1 shows that, for R = 0, unphysical con- 
tributions dominate the F (= N) cross section over most 
of the angular range, i.e. for 9 > 50°. In particular, 
the F curve is completely different from the expected ex- 
ponential decrease. Also the F (= — N) LAM « for 
9 > 50°. At forward angles, oscillations in the F LAM 
curve indicate that another F subamplitude is present, 
which interferes with the unphysical one. This behavior 
does not support the conjecture that the F LAM of this 
other subamplitude has the expected value, A = 10.0. 

By comparing the lower panel of Fig. 1 with the cor- 
responding Fig. 1 of [^, for which A = 10 and A = 2.0, 
one observes that the change in A has not altered the 
magnitude of the unphysical contribution in the angu- 
lar region where it is dominant. The explanation for 
this observation is straightforward. The amplitude ob- 
tained on dropping Si from the term 5*; — 1 in (H) is 
fs{9) = iS{l — x)/k. Applying to f5{9) the same proce- 
dure used by Fuller to derive the NF Coulomb subampli- 
tudes (Ref. Eq. 11), we find that the NF components 

of fsi9) are /i^^(6') = ±[27rfc(l - x)]-\ The correspond- 
ing NF cross sections, downward shifted by one vertical 
unit, are shown in Fig. 1 by the thin dotted curve and la- 
belled NFs- This curve shows that the major part of the 
unphysical contribution, when it dominates the Fuller F 
subamplitude, is due to /j^' (9) and does not depend on 
the properties of Si. 

A more satisfactory Fuller NF result is obtained by 
dropping 1 from the term 5'; — 1 in (^), and then consid- 
ering the PWS thus obtained as a distribution. To ensure 
convergence of the resulting NF subamplitudes, a YRW 
resummation is performed on them. This is explained 
in more detail in Sect. II], and in a systematic nota- 



tion developed there, is denoted R = Oy. Figure 1 shows 
that the R = Oy results are rather good at forward an- 
gles, 9 < 70°. Apart from a small region around 9 = 0°, 
where the condition Asin^ ^ 1 is not satisfied, the N 
(= — F) LAM agrees closely with the expected value of 
—A up to « 70°, and the N (= F) cross section curve 
follows the expected exponential decrease. For 9 > 120°, 
the N cross section is still dominated by an unphysical 
contribution. At intermediate angles, 70° ^ 9 < 120°, 
interference oscillations appear both in the N cross sec- 
tion and in the N LAM curve. It is interesting to note 
that the NF LAMs are more sensitive to interference ef- 
fects in Fig. 1 than are the NF cross sections. Also, in 
the interference region, one cannot attach the meaning 
of a local angular momentum to the subamplitude phase 
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FIG. 1: Strong absorption model N (thick continuous curves) 
and F (thick dashed curves) cross sections (lower panel) and 
LAM (upper panel) calculated using the R = and R = Oy 
NF methods. The thin curve shows the full cross section. The 
thin dotted curve (NF^) shows the F (= N) cross section for 
the unphysical amplitudes fs^\9) (displaced downward by 
one unit). The thin dot-dashed line shows the slope of the 
expected exponential behavior for the NF cross sections 



derivative. In our case, in this interference region, the N 
LAM curve oscillates around the expected semiclassical 
value of —A in the region, 70° ^ 9 < 90°, where the true 
semiclassical component dominates the N subamplitude, 
and around the unphysical value of at larger angles. 



C. Optical model for ^''O 



elastic scattering 



Figure 2 shows our results for the phenomenological 
(WS2) optical potential used to fit |§ the ^^O + ^^O 
elastic cross section at i^iab = 145 MeV. The usual Fuller 
NF method has been applied that employs an analytic 
formula for the NF subamplitudes of the Coulomb scat- 
tering amplitude The parameters for this potential 
are given in Table 1 of Ref. . In the upper panel we dis- 
play LAM/fc, which we call the Local Impact Parameter 
(LIP), and in the lower panel a log plot of a-{9) sm9. The 
thin continuous lines show the (unsymmetrized) cross 
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E,„^=145MeV 



60 120 180 

e(deg) 

FIG. 2: Fuller N (thick continuous curves) and F (thick 
dashed curves) cross sections (lower panel) and LIP (upper 
panel) for ^*'0 + ^"O clastic scattering at Ei^h = 145 McV. 
The thin continuous curves show the cross section and LIP 
using the full amplitude. The two thin dotted lines show 
interpolations of the average behavior of the N cross section 
in the ranges 50° < 6* < 90° and 6 > 150°. The inset in the 
lower panel show a vertical magnification by a factor of 3 of 
the thin rectangular area. 



section and LIP for the full amplitude. For 9 « 20°, 
and at backward angles, the full cross section is charac- 
terized by oscillations having approximately the same pe- 
riod. In the intermediate angular range, 50° ^ < 150°, 
there are oscillations of longer period, which increases as 
6 increases. These properties, which axe typical of light 
heavy-ion scattering, are much more complicated than 
those for the simple SAM cross sections, and indicate 
the presence of several interfering subamplitudes. 

The oscillations in the LIP curve for the full amplitude 
qualitatively follow those in the full cross section. In ad- 
dition, the change in sign of the full LIP for 6 « 20° indi- 
cates a crossover between one N subamplitude, dominant 
at smaller angles (where the full LIP oscillates around 
negative values), with a F subamplitude that is domi- 
nant at larger angles (where the full LIP oscillates around 
positive values) . The oscillations of longer period in the 
intermediate angular range for the full cross section corre- 



spond to an oscillating positive LIP. This indicates there 
is an interference between two F subamplitudes with a 
crossover near the deep minimum at Ri 70°. The oscil- 
lations at backward angles in the full cross section indi- 
cate an interference of one of these F subamplitudes with 
its continuation to the opposite side of the target (N). 

The results of the Fuller NF method (thick curves, con- 
tinuous for N and dashed for F quantities) in Fig. 2 
confirm this interpretation. The N cross section crosses 
the F one at « 20°; the F cross section oscillates in 
the intermediate angular range before decreasing, almost 
monotonically, towards backward angles, where it meets 
the N cross section. For 9 > 60°, the behavior of the 
N cross section is characterized by rather complicated 
oscillations. It is difficult to imagine that these oscilla- 
tions arise only from interference between a N subam- 
plitude, which is dominant at forward angles, with a N 
subamplitude dominant at backward angles. The thin 
dotted lines in Fig. 2 show interpolations (linear fits) of 
the average behavior of the N cross section in the ranges 
50° < 6* < 90° and 9 > 150°. In the crossing region of 
the dotted lines, the N cross section is about two orders 
of magnitude larger then the crossing value. The contri- 
bution from an additional N subamplitude is apparently 
necessary to explain this behavior of the N cross section. 
The shape of the N LIP curve supports this conjecture. 
The N LIP curve is oscillatory, with increasing oscilla- 
tion amplitude, up to ^ ~ 85°, where the deep minimum 
corresponds to a crossover between two interfering N sub- 
amplitudes. For 85° < 6* < 145°, the N LIP curve is oscil- 
latory around 0, which suggests an unphysical origin for 
the N subamplitude dominant in this angular range. At 
9 ~ 145°, where there is a narrow maximum in the N LIP 
curve, the unphysical N subamplitude crosses a different 
N subamplitude, that becomes dominant at backward an- 
gles. Oscillations, similar to those observed at backward 
angles in the N LIP curve, are present at backward an- 
gles in the F LIP curve. The corresponding oscillations 
in the F cross section are barely visible with the scale 
used in the lower panel of Fig. 2. They can, however, be 
observed in the inset, where the thin rectangular area is 
plotted with a vertical scale magnified by a factor of 3. 

In summary our detailed analysis of the NF LIP and 
NF cross sections indicates the presence of an unphysical 
contribution, which makes the properties of these quan- 
tities unneccessarily complicated. 



D. Additional optical model examples 

In this section we briefly discuss flve additional ex- 
amples where the NF optical potential cross sections are 
dominated by unphysical contributions over wide angular 

ranges. 

Figure 3 shows N (upper panel) and F (middle panel) 
cross sections obtained by the usual Fuller NF method 

for different optical potentials that describe the elastic 
scattering of ^'^O by ^^O, at Ei^^ = 480, 704 and 1120 
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titles, logiQ[a{0) sm6], to the values of the corresponding 
N quantities, indicated by N{0q), at 6*0 = 150°. Except 
for a restricted region at forward angles, whose width 
depends on the system considered, the normalized Fuller 
NF cross sections are practically equal, for different col- 
liding partners and/or energies. In every case, the un- 
physical effects of the Fuller NF meth od a re as striking 
as those for the SAM discussed in Sect. II B. In the lower 



60 



120 



9 (deg) 



panel we have plotted, without normalization factors, the 
fuU cross sections and, for 9 > 120°, the Fuller N F) 
cross section. It is very apparent that, at large angles, the 
behavior of the full cross sections is quite different from 
that of the Fuller NF cross sections. In particular, the 
160 + 160 full cross sections at 480 MeV and 704 MeV 
possess a pronounced interference pattern in an angular 
region where the Fuller NF cross sections are dominated 
by unphysical contributions. For completeness, we re- 
mark that the NF LIP is practically null in the angular 
ranges where the NF cross sections show unphysical be- 
havior for all these cases. 



III. IMPROVED NF METHOD USING 
RESUMMATION 

A. Resummation of full and NF PWS 



In the preceding Section, we have discussed examples 
where the Fuller NF method resulted in unphysical con- 
tributions. We also clearly demonstrated that the Fuller 
method has the capability to identify the unphysical con- 
tributions through anomalous behavior in the NF LAM, 
or LIP. This capability helps us to avoid misleading inter- 
pretations of the full and NF cross sections obtained from 
the Fuller NF method. However the problem of finding 
more satisfactory NF methods without unphysical con- 
tributions remains open. 

A possible solution to this problem was proposed by 
Hatchell [|l2|, who used a modified NF method. The 
180 modifications consisted of, first, in writing f(6) in the 
resumnied form (x ^ 1) 



FIG. 3: Fuller N (upper panel) and F (middle panel) optical 
potential cross sections for ^^0 + ^"^O, a -I- ^^C and a + ''"Ca 
at different energies. All the N and F curves are normalized to 
the values of the corresponding N quantities ai 6 = 150° . The 
lower panel shows the full cross sections and, for 9 > 120°, 
the unphysical behavior of the Fuller N (« F) cross sections. 



MeV (potentials labelled WS2 in Ref. §), and of a par- 
ticles by 12 C and 4°Ca, at Ei^^ = 1370 MeV (potentials 
labelled WSa in Ref. [0). We used relativistically cor- 
rected kinematics in all these optical potential calcu- 
lations, and symmetrisation effects were ignored for the 
160 + 160 system. 

In the N and F panels, we have normalized the quan- 



1 1 °° 

^ ' 1=0 

r — 1,2,..., and, second, in using a different splitting 
in the resummed PWS (^) for the Legendre polynomials 
into traveling waves. 

The use of the resummed form ^ for f{0) was orig- 
inally proposed |Q by Yennie, RavenhalL and Wilson 
(YRW) in a different context. Equation ^ is an exact 
resummation formula, of order r, which is derived from 
the recurrence relation for Legendre polynomials. Some 
mathematical properties of the resummed PWS (^ have 
been investigated by Wimp |Q. The YRW resummed 
form (i) can be derived by iterating r times, starting 
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from = a;, the resummation identity [x ^ \) 

oo oo 



i=0 



where 



(0 



I 



21- 1 



h-i 



2Z + 3"'+i ' 



(10) 



(11) 



with a^^^ 



0. 



It is important to remark that the resummed coeffi- 
cients a^j^\ being hnear combinations of the af^ = ai, 
can have very different properties from the original phys- 
ical ai. However, no information about the physical a; 
is lost on applying the resummation procedure, and the 
value of f{9) does not depend in any way on the resum- 
mation order used. This is true for (^) and for all the 
resummed forms for f{9) derived in this paper and is 
a consequence of the (exact) mathematical properties of 
the Legendre polymonials. 

Equation (|^) does not hold for a PWS written in terms 
of a linear combination of Legendre functions of integer 
degree of the first and second kinds (x ^ ±1) 



where 



Ci{x) = pPiix) + qQi{x), 



(12) 



(13) 



with p and q real or complex constants independent of 
I. Rather, as a result of the property lQi-i{x) ^ 1 as 
I —^ 0, the exiendec? resummation identity holds (x ^ ±1) 
PI 



OO ^ 

E«'-^^A(x) = ^ 

!=0 



Li=0 



(14) 



where the recurrence relation ( |l 1| ) gives cif' in terms of 

a|* . By iterating ( p^ ) r times , T(&) can be written in 
the extended resummed form {x ^ ±1) ||25|| 



1=0 



(1 - xY 



(15) 



r of the resummation. In particular, the last term on 
the r.h.s. of (^5|) is omitted from these r dependent re- 
summed NF subamplitudes. However, the sum of the NF 
resummed subamplitudes remains f{0), because the dif- 
ferences exactly cancel each other. A different extended 
resummation identity occurs in the Hatchell approach, 
because the functions he used in place of the Qi (x) obey a 
(inhomogeneous) recurrence relation, different from that 
for Qi{x) (see Ref. for details). 

It is important to realize that the YRW resummation 
procedure (which has been extended in to treat PWS 
like (|l^) ) can also be used to speed up the convergence of 
the PWS (0) or (^), or even to ensure their convergence if 
the PWS was originally defined only in a distributional 
sense: in fact if ai - 0{l-P) then a'f^ - ©(r^-^r) 
Indeed the YRW resummation procedure was originally 
introduced |Q| to produce, for r > 1, a convergent 
p6[ PWS for the full amplitude of high-energy electron- 
nucleus scattering in which Si is asymptotically Coulom- 
bic, and also to speed up the convergence of this PWS. 

Using his splitting for Pjjx) into traveling angular 
waves, Hatchell has shown ||l2[ that the unphysical con- 
tributions to the SAM NF cross sections systematically 
decrease on increasing r, i.e. on using the resummed 
PWS (^ before the splitting of Pi{x), rather then the 
original unresummed PWS (^. More recently [|l6[ pT| , 
it was shown that the same NF resummed method gives 
even better results if the Fuller Q'f^\x) functions are 
used. The superiority of the Qp^x) seems to be con- 
nected with the greater rapidity with which the qI^\x) 
approach their asymptotic behavior (0) |l^, |l^ , com- 
pared to the Hatchell NF functions. 

The success of using (^ before applying the NF split- 
ting (0) depends upon the properties of the resummed 



coefficients a'f^ For the SAM at low / values, the 
rapidly decrease in magnitude [ p^ with increasing r. As 
a result, the low / values, where the splitting into running 
angular waves is physically less reasonable, give a smaller 
contribution to the resummed PWS. 

However, in some cases, this resummation technique 
acts in the opposite direction, by enhancing the undesired 
unphysical contributions to the NF resummed subampli- 
tudes. We have found that this happens, for example, 
for pure Coulomb scattering, for scattering by an impen- 
etrable sphere, and for the SAM (see |]l9j) when the cross 
section is calculated at an angle tt — 0, using the identity 
P; [cos(7r - 61)] = (-l)'P;(cos6'). 



(r) 



It is important to note that J-{9) is independent of the 
value used for r. 

The qJ'^^(x) used to spht Pi{x) in the Fuller NF 
method are a particular case (with p = 1/2, and q — 
±i /tt) of the more general linear combination Ci (x) 

When the splitting (||) is made in the initial (^ or 
resummed PWS (||), different NF subamplitudes are ob- 
tained, i.e. the NF subamplitudes depend on the order 



B. Improved resummation of full and NF PWS 



One possible solution to the intriguing puzzle discussed 
in Section HI A is to regard ( p^ ) as a particular case of 
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the modified resummation identity pi 



oo -.00 

^ar)p,(.) -^^j:4'Pii-i ^ = 1,2,. 



/=0 



(16) 



with a.i + PiX ^ and 

o^' = /^^^-ti'^ + + (17) 

For ai,Pi ^ 0, apart a renormahzation factor, the r.h.s. 
of (^6|) depends only on the ratio Pi/ai. Thus, without 
loss of generality, we can assume ai = 1 for all i. By 
iterating ( p^ ) r times, we can write f(0) in the modified 
resumnied form {1 + Pi ^ 0) 



/ r \ 00 



(18) 



r = 1, 2, . . . . It is straightforward to show that the ex- 
tended modified resunimed form for the PWS (|l^) is 
given by (x ^ ±1 and 1 + Pi ^ 0) 



^ r ^ \ 00 

VJ=1 / /=0 



(19) 



The identity (|l8|) is the key result for our improved NF 
method. The YRW resummation formula (^) is obtained 
when Pi = P2 = . . . = Pr = —1. 

The resummed form ( p^ ) with P = Pi = ... = Pr 
is a particular case of a more general one |l^, which 
uses a basis set of reduced rotation matrix elements; this 
gives the amplitude for more general scattering processes 
than those described by (|]). For these general PWS, a 
Fuller-like NF decomposition can be introduced |2^, 
^ which allows the scattering amplitude to be separated 
into NF subamplitudes. In some cases, the NF cross 
sections contained unexpected (unphysical) oscillations 
p^ , which are enhanced if the generalization of (^) is 
used, but which disappear for an appropriate choice of 
the /^-parameter in the generalization of (p^). 

The considerable successes achieved by the original 
Fuller NF method suggest that the modified resummed 
form (^8|) be used to diminuish unphysical contamina- 
tions to the NF subamplitudes when they are present. 
To do this, we must give a practical rule to fix the values 
of the Pi parameters. In Refs. |l^, it was proposed to 
select the value of /3 = /3i = . . . = /3r so that (1 -I- Px)~^ 
approximately mimics the shape of the angular distribu- 
tion. The shape of the cross section can however be very 
different from that given by (1 -|- Px)~^ . It is therefore 
desirable to test a different recipe, based on a simple rule. 
The quantitative recipe proposed here is inspired by the 
observation that the modified resummation procedures 



produce a more physical NF understanding by reducing 
the contribution from the low I values in the resummed 
PWS. This suggests that we select the /3i, /?2, . . . , /3r in r 
repeated applications of (|l6|), so as to eliminate as many 
low I terms as possible from the resummed PWS in (|T^). 

} to {al''} is linear tridi- 



The transformation from {( 
agonal, with coefficients linear in Pi, which means that 
a resummation of order r allows one to equate to zero 
the leading r coefRcients a] , with Z = 0,1,... ,r— l,by 
solving a system of r equations of degree r in the parame- 
ters Pi,P2, ■ ■ ■ , Pr- We will call the resummation defined 
in this way an improved resummation of order r. 

It is straightforward to show that the improved resum- 
mation of order r = 1 is obtained by choosing 



Pi = -3ao/ai, 
while the improved resummation of order r 

by 



(20) 
2 is given 



f3i 



{B±VB^-4A)/2, 



with A and B solutions of the linear equations 



^aiB — — (2q 
(ao + fa2)B = -ai. 



(21) 



(22) 



Higher order improved resummations require the solution 
of more complicated systems of equations. 

Note that the improved resummation of order r = 1 
is, obviously, not defined if ai — 0. Similarly the resum- 
mation of order r = 2 is not defined for an (accidentally) 
zero value of the determinant of the linear equations (|2^). 
Analogous limitations are expected to hold for higher or- 
der resummations. However, these pathological situa- 
tions were never observed in our calculations. In all the 
cases we have analyzed using r < 2, we find that the 
improved resummations considerably reduce the width 
of the angular regions in which the Fuller NF cross sec- 
tions exhibit unphysical behavior. In these analyses, we 
have used Si from simple parametrizations, as well as 
from some of the optical potentials currently employed 
to describe light heavy-ion scattering. 

Note that, as for the YRW resummed form (^, no 
physical information is lost on using our improved re- 
summed form (|l^), and the full f{9) does not depend 
in any way on the resummation order used. This might 
seem surprising, because our method omits the contribu- 
tion from low values of I. For example, for the improved 
resummation of order r — 1, the resummed term Oq^'' is 
set to zero in the resummed PWS, but the information 
on the physical value of ao is contained in the resummed 
term a^i^ and in the resummation parameter /3i, and sim- 
ilarily for higher order resummations. The improved re- 
summation method only modifies the NF (resummed) 
subamplitudes, not their sum f{9), as it tries to elim- 
inate unphysical contributions that have their origin in 
the NF splitting of the Pi{x) into running waves q['^\x). 
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The physical meaning of the NF spHtting (a) is based 

on the asymptotic properties of Pi{x) and Qf^\x), see 
(0) . The NF sphtting (||) is not expected to be physically 
meaningful for low / values. For example, for Z = 0, (|^) 
states that Po{x) = 1 is the sum of 



g[,^\:r)^i±f ln[(l 



(23) 



This decomposition of unity is undoubtly mathematically 
exact, but it is difficult to think that it can be physi- 
cally meaningful. It is contributions of this type that the 
improved resummation omits from the resummed PWS, 
after having moved to higher values of I in the resummed 
PWS, the physical information contained in the original 
Si. 

It is also interesting to understand the reasons for the 
successes and failures of the Hatchell method for the SAM 

f, |l^ based on the use of the YRW resummed form 
for f(9) and the Fuller NF splitting (||). In those 
cases where the Hatchell method is successfuU, we find 

(r) 

for low I values. Si « Si+i, and the aj are made small in 
magnitude with the choice /3i « /32 ~ . . . ~ /?r ~ — 1- In 
those cases where the Hatchell method fails, it happens 
that, for low / values, Si ~ — and the choice /3i « 

(r) 

/32 « . . . « /3r ~ 1 makes the aj small in magnitude. 
Our improved method automatically chooses appropriate 
values for the resummation parameters using the exact 
values of Si . 

Finally, we summarize the procedure used by our new 
improved resummation method. The calculations are 
performed applying: first, an improved resummation of 
order r — 0,1,2, with r — meaning no resummation, 
second, the Fuller N F sp litting (||), and third an extended 
YRW resummation ( [l5|) of the NF (resummed) subampli- 
tudes. This latter YRW resummation is necessary to en- 
sure the convergence of the final NF PWS. This final re- 
summation can, however, be replaced with any other pro- 
cedure that provides convergence of the final NF PWS. 
The results obtained from these three steps will be in- 
dicated by the notation R — Oy, ly, 2y. The notation 
R = Oy is used in Fig. 1. Also, we use R = to indicate 
the original Fuller method for the SAM case: no resum- 
mation and no final YRW resummation, because, in this 
case, the presence of the 1 in the term Si — I ensures 
rapid convergence of the PWS. 



C. Strong absorption model for elastic scattering 

Our results obtained after ap plyin g the R — ly,2y 
procedures to the SAM of Sec. II B are shown in Fig. 
4. The effectiveness of the improved resummation pro- 
cedure is impressive, both for the NF cross sections and 
for the NF LAM curves. Using the R = ly method 
(for which /?! = -0.7621), the F (= -N) LAM and 
the F {— N) cross section are in agreement with the ex- 
pected results up to « 140°. For R = 2y (which has 
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FIG. 4: Same as Fig. 1, except using the R — ly, 2y resum- 
mation methods. 



/3i,2 — — 0.8634 ± 0.0723 i), the agreement covers almost 
the whole angular range. The small irregular oscillations 
appearing at large 6 for the N LAM curve, when R = 2y, 
may arise from precision limitations (64 bit floating point 
representation) in the calculations, or from residual un- 
physical contributions not completely removed by our im- 
proved method. 

If instead we apply the YRW resummation procedure, 
we find that it gives worse results compared to our R = 
ly, 2y methods, although better results than the Fuller 
R = Oy method. This is because the YRW choice of 
f3i = (32 — —1 is not too far from our improved estimate 
of these parameters. 

Note that, using our improved procedure, the value of 
the SAM cut-off parameter, A = 10, is correctly identified 
in the plots of the LAM. Also, no shift of the LAM oc- 
curs on changing the resummation order. We recall that 
resummation does not change or translate in any way the 
physical properties Si ; nor does it alter the physical con- 
tent of f{9)- It only changes the coefficients aj'^'' in the 



resummed PWS (|18|). Note that the a] have a different 



meaning from the partial wave amplitudes a;, and the 
summation index of the resummed PWS (^8|) should not 
be identified as the orbital angular momentum quantum 
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FIG. 5: Same as Fig. 2 except using the R = 1y NF 
resummation method and the full LIP curve is not shown. 
Cross sections and LIP calculated using classical mechanics 
are shown by thin curves (continuous for N and dashed for F 
contributions). The indices 1, 2, 3 and 4 identify the curves 
corresponding to different branches of the classical deflection 
function. 



number. This is because the resummation extracts the 
factor Jll=i(l + Pi^)^^ from the Legendre polynomial 
PWS. 



D. Optical model for + ^^O elastic scattering 

Figure 5 shows our results for the WS2 optical po- 
tential of Sec. p^I C| using the R = ly method (for which 
/3i = -0.9997- 0.0800 i). The oscillations of the Fuller N 
curves around 6 ~ 90° are removed, and both the R = ly 
N cross section and the N LIP curve are considerably sim- 
pler than those from the Fuller method. The F curves 
are essentially the same as those in Fig. 2, although the 
R = ly method has suppressed the oscillations in the F 
cross section at backward angles (compare the insets in 
the lower panels of Figs. 2 and 5), as well as the oscilla- 
tions in the F LIP curve in the same angular range. This 
demonstrates that the R = ly method is more effective 
than the usual Fuller one in decomposing f{9) into slowly 



varying NF subamplitudes. It also shows that some os- 
cillations in the Fuller NF cross sections in Fig. 2 are 
artifacts, without any physical meaning, introduced by 
the properties of the NF technique used. We have also 
applied the improved resummation method R = 2y. The 
results are practically the same as those for R = ly and 
are not shown. 

The cleaning by the R = ly procedure in Fig. 5 of the 
original Fuller NF subamplitudes (Fig. 2) is impressive 
and allows a better identification, in the F cross section 
and F LIP of the dominance, for 9 > 90°, of the con- 
tributions from classical-like trajectories refracted from 
the internal part of the nuclear interaction. In Fig. 5, 
this interpretation is demonstrated by the overall agree- 
ment, for 6 > 90°, between the F curves and the classical 
mechanical results (thin curve, labelled 1) corresponding 
to impact parameters smaller than the classical orbiting 
one. 

For the WS2 potential the collision energy is, in fact, 
slightly below the critical energy at which orbiting dis- 
appears for the classical deflection function (when it is 
transformed into a nuclear rainbow minimum). Because 
of this, the dependence of the impact parameter b on 
is that of an infinitly many- valued function ([3C[|, p. 
127-129). Four of these branches, corresponding to the 
deflection function being larger than —360°, are plotted 
for the LIP in the upper panel of Fig. 5 (shown by thin 
curves) and labelled from 1 to 4 for increasing values of b. 
The LIP is assumed to be equal to b for the F branches 
(thin dashed curves, with labels 1 and 4) and to —b in 
the N case (thin continuous curves, 2 and 3). In the 
lower panel, the thin curves, with the same labels, show 
the contributions to the classical cross section from these 
branches, in which we have included, in the usual simple 
way (fltl, p. 49), the absorptive effects of the imaginary 
part of the optical potential. 

The agreement between the quantum R = ly F curves 
and the classical curve labelled 1 is overall good for 
6 > 90°, and is impressive for 9 > 120°. This shows 
that the R = ly F subamplitude, calculated by an exact 
quantum method, is dominated for 9 > 90° by a classi- 
cal contribution corresponding to trajectories with small 
impact parameters, 6 < 5 fm, which are refracted by the 
nuclear interaction. The contribution from F branch 1 
continues at 180° into the contribution from N branch 2. 
However, as 9 decreases from 180°, there is an increasing 
disagreement with the average quantum N contribution. 
Evidently, for impact parameters that approach the or- 
biting one (fcorb = 6.001 fm), diffractive effects start to 
become significant. Finally, we note that there is a large 
disagreement between the classical F curve 4 and the av- 
erage behavior of the quantum F cross section at forward 
angles. This suggests that the other F subamplitude re- 
sponsible for the interference pattern in the F quantum 
cross section cannot be considered a classical-like refrac- 
tive contribution. 
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FIG. 6: Same as for the upper and middle panels of Fig. 
3, except without normalization factors and the R = 1y NF 
resummation method has been used. The inset shows the full 
cross sections. 



E. Additional optical model examples 

In Fig. 6, we show the N (upper panel) and F (lower 
panel) cross sections calculated using the improved NF 
met hod, R = ly, for the optical potentials considered in 
Sec. pi D| . In contrast to Fig. 3, no normalization factors 
are used in Fig. 6. Comparison of Fig. 6 with Fig. 3 
shows that the R = ly NF resummation method has 
removed the unphysical contributions that dominate the 
Fuller NF cross sections over most of the angular range. 

For these cases, /3i = -1.0001 - 0.0283^,-1.0001 - 
0.0218^,-0.9991 - 0.0102^,-0.9958 - 0.0021i, and 
—0.9985 — O.OOlOi, in descending order for the collision 
systems listed in the legend of the upper panel of Fig. 
6. In all cases, /3i « —1, which is a consequence of the 
fact that, for these five collision systems, at low I values 
we have Si ~ Si+i. The result /3i « — 1 also explains 
why the unphysical contributions in the Fuller NF cross 
sections of Fig. 3 have very similar shapes. Indeed for 
/3i = —1, the angular dependence of the unphysical cross 
sections is oc (1 — x)~^ oc sin"'' (0/2), which is the same 
as the Rutherford cross section. This implies that on 
making the (usual) plot of the ratio of the cross sections 



to the Rutherford one, the unphysical contributions will 
appear independent of 0, i.e. a constant. 

Figure 6 clearly shows that the R = ly method has 
cleaned the Fuller NF cross sections of unphysical contri- 
butions. The R = ly and Fuller NF cross sections agree 
closely at forward angles, with the R = ly curves pro- 
viding the correct continuation of the Fuller NF curves 
to larger angles, where the unphysical Fuller NF contri- 
butions become dominant over the R = ly results. In 
addition, the R = ly procedure clearly identifies the os- 
cillatory pattern of the '^O -I- ^^O full cross sections at 
Eiah — 480 MeV and 704 MeV as an interference between 
the NF resummed subamplitudes. 

As a final curiosity, we note that a 'suspect' behavior 
appears in the R = ly NF cross sections for a-f^^C at 
e > 40°. We find that the corresponding R = ly NF 
LIP are « for 6* > 60°. However, the LIP for the full 
amplitude is also « for 9 > 60°, indicating that the 
full cross section is dominated by contributions from low 
partial waves in this angular region. 

We have repeated the calculations, substituting the 
WS form factors used in [0, with symmetrized WS-like 
form factors defined by 



fsym{r, R, d) 



sinh — /(cosh — 
a a 



cosh ■ 



(24) 



where R and d arc the radius and diffuseness parameters, 
respectively. We find that the full cross section decreases 
for 6 > 60° by more than 5 orders of magnitude, which 
exceeds the precision limits of our optical potential com- 
puter code. At forward angles the effect of the substitu- 
tion is very small. This supports the conjecture that the 
suspect behavior in the R = ly NF and full cross sec- 
tions arises from the 'cusp', which the usual WS poten- 
tials have at the origin. The cusp is expected to produce 
diffraction scattering for the low partial waves, which is 
equally distributed between the N and F subamplitudes. 
It is this effect that we observe, and which disappears 
on removing the cusp. Similar effects, although masked 
by the limited range of the scale along the ordinate, are 
also observed in the '^O -I- '^O collision at E'lab — 1120 
MeV, as well as tentatively in the R = ly N a+'^'^Ca 
cross section. 



IV. CONCLUSIONS 

Our new NF resummation procedure clearly improves 
the original Fuller NF method, as is evident from the 
seven examples discussed here, and from those presented 
elsewhere In all these cases, we obtain NF re- 

summed cross sections that are more slowly varying and 
less structured than the Fuller ones. On the one hand, 
our results confirm the utility of NF methods for gaining 
insight into the properties of the subamplitudes respon- 
sible for complicated structures in cross sections. On the 
other hand, they remind us of the empirical origin of NF 
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methods, and suggest eaution in the interpretation of re- 
sults obtained from NF techniques. Using different NF 
methods lets us check what parts of the resulting NF sub- 
amplitudes are independent of the particular technique 
used. Only properties stable with respect to different 
NF methods, can be considered as manifestations of some 
physical phenomenon. 

In addition, we have shown that in all NF analyses, it 
is desirable to investigate the behavior of the LAM. This 
quantity is more sensitive to interference effects than are 
the NF cross sections, and a null value (or an oscillatory 
behavior around zero) of the NF LAM in a certain angu- 



lar range may indicate the dominance of an unphysical 
contribution. 
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